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We use Big Bang Nucleosynthesis calculations and light element abundance data to constrain 
the relative variation of the deuteron binding energy since the universe was a few minutes old, 
8Q — Q(BBN) — Q(present). Two approaches are used, first treating the baryon to photon ratio, 
rj, as a free parameter, but with the additional freedom of varying 5Q, and second using the WMAP 
value of n and solving only for SQ. Including varying Q yields a better fit to the observational data 
than imposing the present day value, rectifying the discrepancy between the 4 He abundance and 
the deuterium and 7 Li abundances, and yields good agreement with the independently determined 
tjwmap- The minimal deviation consistent with the data is significant at about the 4-<r level; 
SQ/Q = —0.019 ± 0.005. If the primordial 4 He abundance lies towards the low end of values in the 
literature, this deviation is even larger and more statistically significant. Taking the light element 
abundance data at face- value, our result may be interpreted as variation of the dimensionless ratio 
X = tUs/Aqcd of the strange quark mass and strong scale: SX/X — (1.1 ± 0.3) x 10 -3 . These 
results provide a strong motivation for a more thorough exploration of the potential systematic 
errors in the light element abundance data. 
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I. INTRODUCTION 

Recent astronomical data suggest a possible variation 
of the fine structure constant a — e 2 /Tic at the 10~ 5 level 
over a time-scale of 10 billion years, see pj (a discussion 
of other limits can be found in Ref. Q and references 
therein). Naturally, these data motivated more general 
discussions of possible variations of other constants. Un- 
like for the electroweak forces for the strong interaction, 
there is generally no direct relation between the coupling 
constants and observable quantitites. In recent papers 
dHHl, we presented general discussions on the possi- 
ble influence of the strong scale variation on primordial 
Big Bang Nucleosynthesis (BBN) yields, the Oklo natu- 
ral nuclear reactor, quasar absorption spectra and atomic 
clocks. Here we continue this work, concentrating on 
BBN. 

One can only measure variations of dimensionless pa- 
rameters. Big Bang Nucleosynthes is sensitive to a 
number of fundamental dimensionless parameters includ- 
ing the fine structure constant a , Aqcd /M-Piank and 
m q / Kqcd where m q is the quark mass and Aqcd is the 
strong scale determined by a position of the pole in the 
perturbative QCD runnung coupling constant. In this 
work we search for any possible variation of m q j Aqcd 
because there is a mechanism which provides a very 
strong sensitivity of BBN to this parameter. 

The first and most crucial step in BBN is the process 
p + n — > d + 7. The synthesis starts at t > 3 sec. when 
the temperature goes down below T < 0.6 MeV and lasts 
until t < 6 min. when the temperature becomes T < 0.05 
MeV. The reaction rate for the above process defines all 
subsequent processes and final primordial abundances of 



light elements. Amongst the factors that can influence 
the reaction rate, the most significant seems to be a vari- 
ation of the deuteron binding ener gy ( this variation was 
discussed in Refs. 0, H S H HM EH)- Indeed, the 
equilibrium concentration of deuterons and the inverse 
reaction rate depend exponentially on it. Moreover, the 
deuteron is a shallow bound level. Therefore the rela- 
tive variation of the deuteron binding Q is much larger 
than the relative variation of the strong potential U , i.e. 
SQ/Q >> SU /U. As a result the variations in the strong 
interaction may be most pronounced via the deuteron 
binding energy. We also take into account the effect of 
variation of the virtual level in the neutron-proton sys- 
tem, which is even more sensitive to the variation of the 
strong interaction. 

The question we address here is whether or not ex- 
isting observations of the primordial abundances of the 
light elements suggest any change in the deuteron bind- 
ing energy at the time of BBN. 

To do so, we use a compilation of light element abun- 
dance data from the literature for 4 He, 7 Li and D/H. As 
we show later, the currently greater experimental preci- 
sion on 4 He results in that element dominating our re- 
sults. The other 2 light elements nevertheless provide 
important consistency checks. 

The data we use for 4 He is presented in Table [I] and 
comprised 14 surveys giving estimates for the primor- 
dial value, Y p , derived using, or by extrapolation to, low 
metallicity in each case. There is clear evidence for sig- 
nificant scatter amongst these 14 values, presumably due 
to unquantified systematics, or if not, intrinsic inhomo- 
geneities. The dominance by 4 He, or indeed by any single 
element, unfortunately increases susceptibility to system- 
atic errors, and we have therefore attempted to explore 
the effect of these in several ways. 

Firstly, in order to make best use of all the available 
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TABLE I: Data on the primordial 4 He mass fraction 





Ref. 


0.2391 


± 0.0020 


[12] 


0.2384 


± 0.0025 


[13J 


0.2371 


± 0.0015 


LLL 


0.2443 


± 0.0015 


[15j 


0.2351 


± 0.0022 


[16] 


0.2345 


± 0.0026 


[11] 


0.244 


± 0.002 


[18] 


0.243 


± 0.003 


[19] 


0.232 


± 0.003 


ronl 


0.240 


± 0.005 


[21] 


0.234 


± 0.002 


[22] 


0.244 


± 0.002 


[23] 


0.242 


± 0.009 


[24] 


0.2421± 0.0021 


[25] 



4 He data, we add a constant term to each of the statistical 
errors on Y p , such that that the normalised x 2 for all 14 
points about the weighted mean value is equal to unity. 
This approach is equivalent to the assumption that all 14 
estimates of Y p are unbiased and Gaussian distributed, 
but that there is an additional systematic component to 
the statistical error which is different (and hence random) 
for each estimate. 

Second, as shown later, smaller values of Y p are less 
consistent with SQ/Q = than larger values. Thus we 
carry out a re-analysis using a subset of the Y p 's, taking 
only the highest values such that the normalised \ 2 about 
the weighted mean value is equal to unity, without in- 
creasing the individual errors by a constant, as described 
above. This procedure selects 9 values from the original 
14. In doing this, we are exploring the consequence of 
there being strong systematics for the small Y p 's, and 
little or none for the high values. This is conservative, in 
the sense that we are minimising our estimate for 5Q/Q. 

Finally, in order to obtain some estimate of the plau- 
sible range on our estimate of SQ/Q, we perform the 
converse analysis, subsetting the data by discarding high 
values of Y p , again such that the normalised % 2 about 
the weighted mean value is equal to unity. This leaves 9 
points. The two samples thus overlap. 

The data on deuterium abundances D/H from quasar 
absorption systems were selected according two criteria: 

(i) Metallicity must be low, so as to more closely reflect 
primordial value: [Si/H] or [O/H] less than or equal to 
-2.0. 

(ii) Must be detection, not upper limit. 

These requirements leave only five data points listed in 
Table HU 

The data for Lithium primordial abundance are shown 
in Table ITTT1 Here A = Log(Y Li ) + 12. 

Applying the first procedure described above, in order 
to obtain /N = 1, we have to add to the individual 



TABLE II: Data on the primordial deuterium abundance 



QSO 


z(abs) 


D/Hxl0" 5 


[Si/H] 


Ref. 


Q1009+299 


2.504 


4.0 ± 0.65 


-2.53 


[26] 


PKS 1937- 1009 


3.572 


3.25± 0.3 


-2.26 [O/H] 


[2JJ 


HS0105+1619 


2.536 


2.5 ± 0.25 


-2.0 


[28] 


Q2206-0199 


2.076 


1.65 ± 0.35 


-2.23 


[29] 


Q1243+3047 


2.526 


2.42 +0.35 - 0.25 


-2.77 [O/H] 


w 



TABLE III: Data on the primordial Li/H abundance 



A 


Ref. 


2.09 +0.11-0.12 


[21] 


2.35 ± 0.1 


[32] 


2.36 ± 0.12 


[33] 


2.34 ± 0.056+0.06 


[3J 


2.07 + 0.16 - 0.04 


[35] 


2.22 ± 0.20 


[36] 


2.4 ± 0.2 


[22 


2.5 ± 0.1 


W 



a's 0.0017 for helium points, 0.344 x 10 -5 for deuterium 
points, and 0.028 for lithium points. For the weighted 
mean values we obtain 

Y p = 0.2393 ±0.0011, (1) 

Yd = (2.63 ± 0.31) x 10~ 5 , (2) 

and 

A = 2.315 ± 0.051. (3) 

The latter value corresponds to the following lithium 
abundance 

Y Li = (2.02 ± 0.22) x 10~ 10 . (4) 

The second and the third procedures are meaningful 
only for the helium points. The number of deuterium 
points is too small and the lithium data points are the 
least scattered. We need only 20% increase in individual 
uncertainties to bring x 2 /N to 1 for the lithium data. In 
addition, the deuterium and the lithium data do not pro- 
duce a significant contribution in determination of 5Q/Q 
which is entirely dominated by the helium data due to 
their high accuracy. 

Keeping 9 upper points for the helium mass fraction 
data, that give x 2 /iV = 0.94, we obtain for the weighted 
mean value 

Y p = 0.2424 ±0.0008. (5) 
If we keep 9 lower points, we obtain 

Y p = 0.2363 ±0.0008, (6) 
which is significantly lower than both in Eq.JSJ) and Eq. 
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II. THE BBN EQUATIONS 

We use the standard BBN set of equations that de- 
scribe the time development of the abundances of the 
elements in an expanding Universe |39| 



(7) 




= -3H, 



Pt = -3H(p T +p T ), 



(8) 



(9) 



j,k,i \ 



Y N t Y Nk 



Y N% Y Nj ' 



13 Ni\N k \ ij ^ kl NtlNA 



(10) 
(11) 



where n B is the density of baryons, is the abundance 
of the element Ai Zi . The right-hand side of Eq. (|10|) cor- 
responds to a reaction 

Ni{ A *Zi) + N^tZj) <-> N k ( A *Z k ) + N^'Zt). (12) 

Pt and pr denote total energy density and pressure, re- 
spectively, 



Pt — P 1 + Pe + Pv + Pb, 



PT =P-y + Pe+Pu +PB- 



(13) 



(14) 



Eq.0 defines the expansion rate. Eq.© defines the 
change in time of the baryon density, and the rate equa- 
tion Eq. I |10| l defines the time evolution of the abundances 
and their final values after freeze-out. The last Ea. (|llfl . 
where n_ and n + are the densities of electrons and 
positrons, is the condition of electro-neutrality that de- 
fines a chemical potential of electrons. 



III. EFFECT OF THE DEUTERON BINDING 
ENERGY VARIATION 



change the position of the virtual level e v as well. The 
relation between SQ and Se v can be obtained using the 
fact that both a real level and a virtual one are close to 
E = 0. The relation is (see Appendix) 



SQ 



(15) 



The cross section for n+p — > rf+7 reaction can be found 
in textbooks [ijj . In the leading order in Q /e v the prod- 
uct of the cross section and the velocity is proportional 
to 

gv ~ Q 5/2 /e v . 

Thus, in linear order in SQ we have the following modi- 
fication of the reaction rate 



np — >ci7 




(16) 



We should note, however, that according to our BBN cal- 
culations the direct effect of the deuteron binding energy 
variation (due to the exponential dependence of the in- 
verse reaction rate d'y — > np) is more important than the 
variation of the cross-section np — ► dry Eq. (|16fl . 

We modified one of the standard BBN codes |4l| in 
such a way that Q can be changed for this reaction. 
Varying Q changes the abundances of all three elements 
under discussion. In Fig.l we plot the abundance of D, 
the mass fraction of 4 He, and the abundance of 7 Li, as 
functions of Q at the value of the baryon to photon ra- 
tio 77 = 6.14 x ltU 10 found from anisotropy of cosmic 
microwave background [42| From Fig 1. we see that the 
deuterium abundance is not very sensitive to Q. The 
data are fully compatible with the present value of the 
deuteron binding energy. Such a poor sensitivity can be 
explained by relatively large error bars for the deuterium 
abundance. 

The data on 4 He, in contrast, show strong sensitivity to 
the deuteron binding energy favoring for lower Q during 
primordial nucleosynthesis. The data on 7 Li also favoring 
for lower Q approximately for the same SQ as 4 He. 

The above Figure 1 give a qualitative picture of the 
dependence of light element abundances on the deuteron 
binding energy. In order to obtain more quantitative re- 
sults we analyse the likelihood functions as functions of 
Q and 77. 



The sensitivity of the reaction rate T^d^pn to param- 
eters of the strong interaction in general, and to the 
deuteron binding energy in particular, comes from two 
sources. First, the reaction rate T 1( i^pn depends ex- 
ponentially on the deuteron binding energy Q. Second, 
the cross section of the reaction np — > T-yd is very sensi- 
tive to the position of the virtual level with the energy 
e v = 0.07MeV. Any change in the strong NN-potential 
causing a shift in the deuteron binding energy Q will 



A. The likelihood functions 

The likelihood function for the abundances have been 
choosen in the form 

Lf( v , Q) = cx P (-i J2( Y ^ h ^ Q) ~ F i e3 >« 

x(Y J th (n,Q)-Y j n). (17) 



4 
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QxlO 9 (K) 

FIG. 1: The light elements abundances as a functions of the 
deuteron binding energy Q for rjwMAP = 6.14 x 1(T 10 . The 
vertical line shows the present value of Q = 25.82 x 10 9 K. 
The shaded regions illustrate the la range in the data. For 
helium the high Y p value (EqQ>J is shown. 



Here the sum goes over three light elements, u>ij is the in- 
verse error matrix that was calculated using the approach 
proposed in Ref. 0] . The errors in theoretical values of 
the abundances can be found from the uncertainties in 
the reaction rates 



5Y 4 



Hi 



Y t th J2 A « 



' Rk ' 



(18) 



where ARk are the reaction rate errors, and 
A* 



are the logarithmic derivatives. The error matrix er^ can 
be calculated then by 



j 2 - 



k 



fAR k 



V Rk 



(19) 



The uncertainties in the experimental data 0,(01,0 
should be added to the diagonal matrix elements of the 
error matrix l|19fl 



2 , RrE 2 

cr,-, + tr,- 



(20) 



For 4 He a tot differs from a ex insignificantly, while for D 
and especially for 7 Li an and a\ x are comparable. If we 
neglect the correlations then the matrix ro.y is diagonal 
and equal to 

Wu — 1/CTj 

In this case we can present the likelihood function l|17|) 
as a product of three individual functions Lf(rj, Q) = 
LfD{ r n,Q)Lf H e{r),Q)Lf Li (i-i,Q). The equations 



Y i ( V ,Q) = Y? 



(21) 



defines three lines in 77 — Q plane where the individual 
likelihood functions are equal to one. And the equations 



(Y l ( Vl Q)-Yr) 2 



tot 2 



(22) 



dlnR k 



define la ranges around these lines for each element. 
These ranges are shown in Fig 2. The slope of the deu- 
terium range is smaller than that of helium and lithium 
reflecting smaller sensitivity in Q and higher sensitivity 
in 77. . 

In contrast, the helium range goes almost vertically re- 
flecting high sensitivity of the helium fraction to Q and 
low sensitivity to 77. This low sensitivity to 77 can be ex- 
plained by a large helium binding energy. Only gamma's 
with the energy E 1 >20 MeV can significantly change 
the number of helium nuclei. At any 77 the number of 
such 7-quanta is small at the BBN temperature. We 
can, therefore, expect the low sensitivity of the helium 
mass fraction to 77. 

The Lithium range has two distinct branches cor- 
responding to two different solutions of Ea. (|21() for 77 
at given Q. All three ranges intersect near 77 = 6.5 
and Q = 25. One can expect that the general like- 
lihood function l|17|) will have a maximum in this re- 
gion. Indeed, we found the maximum of Lf(jj, Q) at 
the point 7j m = (6.51 + 0.77 - 0.66) x 10~ 10 and 
Q m = (25.26 ± 0.20) x 10 9 K. Fig. 3 shows la elliptic 
boundary near the maximum. The long axis of the ellip- 
sis is almost vertical. Therefore, the correlation between 
A77 and AQ is not significant. Comparing Fig. 3 and 
Fig. 2 one can conclude that the error AQ is determined 
mostly by 4 He mass fraction data. It is interesting to note 
that 77 m is compatible with the one found from recent 
CMB anisotropy measurement 01 ■ The dark shadow re- 
gion shows the la range for 77 fitted from BBN only at 
present value of Q = 25.82 K. 
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23 24 25 26 27 



Q*1(T (K) 

FIG. 2: lcr-ranges around the maxima of individual likelihood 
functions. The solid lines show lcr-ranges for D, the dashed 
lines are for 4 He (using Y p from EqQ>J, and the dot-dashed 
lines are for 7 Li. For lithium, there are 2 solutions for 77 and 
Q, hence the shape of the error contours is more complicated. 

B. Constraint from CMB anisotropy measurements 

The value of r\ found from CMB anisotropy measure- 
ments 

r)o = (6.14 ±0.25) x 1CT 10 

has rather high accuracy. It is natural to use the 
constraint from this measurement in our study of the 
deuteron binding energy effects. To do this we construct 
another likelihood function which is a function of Q only. 

HQ) = f°° cM- iri 2 J° )2 )Lf(v, Q) dr,- (23) 

If we neglect nondiagonal elements in uiij we can con- 
struct the individual likelihood functions for D, 4 He, and 




24 24.5 25 25.5 26 



QxlO 9 (K) 

FIG. 3: lcr-range about the maximum of Lf(rj, Q) (again us- 
ing Y p from Eq|SJ. The lighter shaded region shows CMB- 
WMAP data for rj. The darker shaded region is the lcr-range 
for t) from BBN calculations using the present-day value of 
the deuteron binding energy, Q — 25.82. A lower value of 
Y p will produce a larger deviation between the t/wmap and 
t/bbn- 



7 Li. They are constructed in the same way as (|23f) using 
instead of general function L/(?7,Q) the individual ones 
LfD{r],Q), Lf He (r),Q) , Lj Ll {r\,Q). These functions are 
plotted in Fig. 4 together with the general likelihood 
function 

From the deuterium likelihood function we found the 
position of the maximum and la deviations: 

Q D = (25.74 + 0.92 - 0.68) x 10 9 . (24) 
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FIG. 4: Individual likelihood functions (12MI for the light el- 
ements. From top to bottom: D, 4 He (Eq|^J, Li, and the 
combined datset. 



The shape near the maximum is apparently non- 
symmetric. The position of the maximum is fully com- 
patible with the present value of Q = 25.82 x 10 9 K. The 
helium likelihood function is much narrower (see the sec- 
ond panel from the top). It gives for the maximum and 



for the la the values 

Q He = (25.37 ±0.13) x 10 9 . 



(25) 



This value lies below the present value of the binding 
energy. Finally, the lithium likelihood function has the 
maximum at 



Q Ll = (24. 



0.43 - 0.59) x 10 9 



(26) 



The position of this maximum is compatible with the 
helium result. 

The general likelihood function (|23|l is plotted in the 
lower panel in Fig. 4 The position of its maximum dif- 
fers only slightly from the position given by the helium 
likelihood function. 



Qbbn = (25.34 ±0.12) x 10 9 



(27) 



It is interesting to compare the light element abun- 
dances for two values of the deuterium binding energies. 
In Fig. 5 we plotted the traditional curves for the light 
element abundances as a functions of r\ for two values 
of Q. The dotted lines in the figures correspond to a 
present value of Q pr esent — 25.82 x 10 9 K, while the solid 
curves correspond to a new value Qbbn = 25.34 x 10 9 
K. Clearly, the new value Qbbn moves the curves closer 
to the data. 

The result which we obtained may be presented as 



SQ/Q = -0.019 ± 0.005, 



(28) 



where SQ = Q bb n — Qpresent ■ If we do not fix rj and try 
to fit it simulaneously with Q we obtain 

SQ/Q = -0.022±0.008, r) = (6.51 + 0.77-0.66) x 10~ 10 . 

(29) 

The obtained rj is fully compatible with the one measured 
by WMAP. 

These values of 5Q/Q and rj were obtained for high 
value of the helium mass fraction Y p . If we use as an 
input the low value of Y p from © we obtain 

SQ/Q = -0.048 ± 0.004, (30) 

If we fit both 8Q/Q and rj we obtain 

SQ/Q = -0.059 + 0.007, r, = (7.55 + 0.91-0.75) x 10~ 10 . 

(31) 

Finally if we use the value of Y p for 4 He obtained using 
the whole sample of 14 points, with increased error bars, 
from Eq.JTJ, we obtain 



SQ/Q = -0.033 ±0.006, 
and for SQ/Q and r\ 



(32) 



SQ/Q = -0.042 + 0.009, rj = (7.00 + 0.85-0.72) x 10" 10 . 

(33) 

The results given in eas. (|28|l and (|30|l therefore represent 
an estimate of the plausible range in SQ/Q. Despite the 
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Combining eqs. il-il'l) and (|34(l we obtain 
S(m s /A Q cD) 



rjx 10" 



FIG. 5: The predicted light element abundance yields as a 
function of rj, for two values of the deuteron binding energy, 
Q. The dotted curve corresponds to the present value of 
Qpresent = 25.82 x 10 9 K. The solid curve corresponds to the 
new value of Q = Qbbn = 25.34 x 10 9 K. The vertical line 
corresponds to rj = 6.14 (WMAP value). The shaded regions 
is the ler-ranges for the observed light element abundances, 
where Y p is from eq|S] 



(tUs/Aqcd) 



(1.1 ±0.3) X 10 



-3 



(35) 



This equation may contain an additional factor (close 
to one) reflecting unknown theoretical uncertainty in 
eq. (|34|l . Note that we obtain here variation at the 
level 10 -3 while the limits on variation of a 0, \^ and 
AgcD/Mpiank 9-13 are an order of magnitude weaker. 
This may serve as a justification of our approach. 



IV. CONCLUSION 

Allowing the deuteron binding energy, Q, to vary in 
BBN appears to provide a better fit to the observational 
light element abundance data. Varying Q simultane- 
ously does two things; it resolves the internal inconsis- 
tency between 4 He and the other light elements, and it 
also results in excellent independent agreement with the 
baryon to photon ratio determined from WMAP. (Fig. 
5). However, the magnitude of the variation is sensi- 
tive primarily to the observed 4 He abundance, which has 
the smallest relative statistical error. A systematic error 
in the abundance of 4 He could imitate the effect of the 
deuteron binding energy variation, although one needs 
a systematic error which is very much greater than has 
been claimed in the most recent observational work. 

We note that Izotov and Thuan |25|. the most recent 
estimate for Y p in our sample, argue that systematics 
are at most 0.6% for that survey. On the other hand, 
the possibility has also been explored that the creation 
of 4 He in population III stars might mean that the true 
primordial 4 He abundance is lower even than that seen in 
the most metal-poor objects . If so, the significance of 
the deviation of SQ/Q from zero we report in this paper 
would be even larger. 

These results hopefully provide an extremely strong 
motivation to obtain substantially better measurements 
of all the light elements, and to explore even more inten- 
sively, the possible sources of systematic errors. 



clear systematic uncertainties in the 4 He data, and ac- 
cepting the WMAP value of rj as being correct, SQ/Q 
appears to deviate from zero by 4cr (eq. I28f) or greater 
(eqs. HUE!!). 

The deuteron binding energy depends on the strong 
scale and quark masses. It is convenient to assume that 
Aqcd is constant , and the quark mass is variable. This 
only means that we measure all energies in units of Aqc d 
(and cross-sections in units Aq^, d ). In Ref. p| we con- 
cluded that the deuteron binding energy is very sensitive 
to variation of the strange quark mass m s |47| : 

S(Q/A QCD ) = l7 5{m s /A QCD ) 
(Q/^qcd) (m s /A QCD ) 
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APPENDIX 

Let Uo(r) be a critical depth potential for which the 
binding energy Q = 0, and Ut{r) a potential for a proton 
neutron system in a triplet state producing a deuteron 
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with small binding energy Q = 2.22 MeV. If we add to 
the deuteron Hamiltonian a perturbation 

5U x {r) = \(U (r) - U t (r)), 

then, variation of A from to 1 will move the binding 
energy Q from 2.22 MeV to 0. From a virial theorem for 
a quantum system we have 



dE_ 



(U (r)-U t (r))x 2 (r)dr, 



(36) 



where x( r ) is the radial s-wave function. For simplicity 
we neglect the d-wave contribution. For Q — ► the main 
contribution into normalization integral for x( r ) comes 
from the region outside of the nuclear forces radius R. 
The normalization integral can be presented as sum of 
contributions from inner and outer regions 



X 2 (r)dr + b 2 e~^ r dr = 1, (37) 
o Jr 



where 7 = y — |j— '-. At \E\ — > the second integral 
dominates giving b 2 = 27. Separating the independence 
of the normalization factor in x( r ) we can rewrite Ea. l36|) 
as 



dE 



d\2 




(f/ (r) - U t (r)) X 2 {r)dr, (38) 



where x( r ) i s practically independent on E inside the 
potential well (where E << U) and x( r ) — e~ 2lT — > 1 at 
r > R when \E\ — *■ . Integrating the left hand side of 
Ea. (|38|l over E from — Q to and the right hand side of 



Ea. (|38() over A from to 1 we obtain 



Q = 



(U (r)-U t (r))x 2 (r)dr 



(39) 



The Ea. (|39|) shows that the position of a shallow bound 
level depends quadratically on the difference between the 
actual depth of the potential and the critical one. For a 
square well Q = Z^^,^ = ^ 

In fact, the Ea. H39|) is valid not only for the energy of 
a bound level but for the energy of a virtual level as well. 
The integration in Eq. (|39|) is over the region r < R where 
the function x 2 (r) is insensitive to the energy E, and the 
quadratic dependence on U — Uq guarantees the validity 
of Eq.® for both U < U and U > U . Thus, for the 
energy of the virtual level we have 



(40) 



TT {U Q {r)-U s {r))r{r)dr 



where U s (r) is the potential for a singlet states. We have 
both Q << U and e v << U. This means that the dif- 
ference between the triplet and singlet potentials is not 
large. Assuming that the changes in the triplet and sin- 
glet potentials are the same we obtain for the changes in 
Q and e v the relation 



Sty 



SQ 

VQ' 



(41) 



This equation also holds for the effect produced by vari- 
ation of the proton mass ( the dominating effect comes 
from variation of Uq). 
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